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Abstract 

Recent research has established that nonsymmetric gravitation theories like 

Moffat's NGT predict that a gravitational field singles out an orthogonal pair 

of polarization states of light that propagate with different phase velocities. 

We show that a much wider class of nonmetric theories encompassed by the 

Xg formalism predict such violations of the Einstein equivalence principle. 

This gravity-induced birefringence of space implies that propagation through 

a gravitational field can alter the polarization of light. We use data from 

polarization measurements of extragalactic sources to constrain birefringence 

induced by the field of the Galaxy. Our new constraint is 10 8 times sharper 

than previous ones. 
04.80,+z, 04.50.+h, 96.60.Tf 
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I. INTRODUCTION 



Gabriel et al. recently established that members of a certain class of nonmetric grav- 
itation theories predict that space is anisotropic and birefringent. These nonsymmetric 
theories, for which Moffat's NGT |J is the prototype, have been studied extensively as po- 
tentially viable alternatives to general relativity, but they predict that a gravitational field 
singles out an orthogonal pair of linear polarization states of light which propagate through 
the field with different phase velocities. The difference in phase velocities could, in prin- 
ciple, be measured in local test experiments [Q] and, so, violates the Einstein equivalence 
principle. In practice, the sharpest constraints on the magnitude of such gravity- induced 
birefringence will be inferred from limits on the cumulative effect that propagation through 
a gravitational field has on the polarization of light. Gabriel || et al. consider the effect 
that the Sun's field could induce and use measurements of the polarization of solar spectral 
lines to impose sharp new constraints on nonsymmetric gravitation theories. Krisher 0] 
considers the effect the Galaxy's field could induce and uses pulsar polarization observations 
to impose a complementary constraint on NGT. 

In this paper we show that the class of nonmetric gravitation theories which predict 
spatial anisotropy or birefringence is far more extensive than the class of nonsymmetric the- 
ories. Indeed, one must take pains to define a nonmetric coupling between the gravitational 
and electromagnetic fields that does not induce anisotropy and birefringence. Consequently, 
searching for effects of propagation through a gravitational field on the polarization of light 
provides tests of the Einstein equivalence principle that have a power, generality and signifi- 
cance comparable to those of more familiar atomic anisotropy (Hughes-Drever) tests || and 
Eotvos and gravitational redshift tests . Recent theoretical developments provide a strong 
motivation for precise new tests of the Einstein equivalence principle. For example, interest 
in scalar-tensor and multi-tensor gravitation theories has been rekindled by the proliferation 
of scalar and tensor fields in effective field theories derived from string physics and several 
threads of quantum gravity research suggest that connection rather than metric may be 
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fundamental. 

The view that searching for evidence of gravity-induced birefringence tests the Einstein 
equivalence principle is not unprecedented. Ten years ago Ni [7] remarked that nonmetric 
gravitation theories encompassed by his \g formalism could predict birefringence of the type 
we study. He also noted that pulsar polarization observations could constrain this possibility. 
Perhaps the significance of these comments was overlooked at the time because no gravitation 
theories predicting such birefringence were known and because the effective time resolution 
of pulsar polarization observations is limited. In any event, their significance is now clear. 
Theories like NGT do predict gravity-induced birefringence and we use galaxy polarization 
observations to constrain this possibility far more stringently than Ni or Krisher were able 
to using pulsar polarization observations. Indeed, the linear and circular polarization data 
for galaxies that we analyze have an effective time resolution that is more than 10 8 times 
that of Krisher's pulsar polarization data. 



II. ELECTRODYNAMICS IN A BACKGROUND GRAVITATIONAL FIELD 

A theory of gravity does two things. It specifies a free, though generally nonlinear, 
dynamics of the gravitational field and it specifies the coupling between the gravitational 
field and matter. The latter determines both how matter generates and how it responds to 
gravity. 

In the case of general relativity these conceptually distinct components of gravitation 
physics are discernible in the standard form of the Einstein field equation, 

G„ v = 8vrT^. (1) 

The structure of the Einstein tensor, on the left-hand side, specifies the general relativistic 
dynamics of the gravitational field itself. The form of the stress-energy tensor, on the right- 
hand side, specifies the way in which matter acts as a source of gravity. 

These same aspects of general relativity and other Lagrangian-based theories of gravity 
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are even more clearly discernible in the structure of the action principles from which their 
field equations follow, 

5 J (C G + C NG ) d A x = 0. (2) 

The gravitational part of the Lagrangian density, Cq, is distinguished by the fact that it 
depends only on gravitational fields. Its variational derivative with respect to these fields 
defines their dynamics. The variational derivative with respect to the metric of general rel- 
ativity's gravitational Lagrangian density, y/—gR/16n, is yf—gG^/lQix. The corresponding 
derivative of its nongravitational Lagrangian density is ^J—gT^jl and, so, the form of Cng, 
which depends on both matter and gravitational fields, specifies the way in which matter 
generates gravity. Its reciprocal role of specifying the way in which matter responds to a 
gravitational field follows from the action principle (0) on variation with respect to matter 
fields rather than gravitational ones. This yields the matter equations of motion. In general, 
these need not follow directly from the gravitational field equations as they do in general 
relativity. 

Experimental tests of the Einstein equivalence principle focus on the behavior of mat- 
ter responding to a background gravitational field. Their significance lies in the way they 
constrain the form of C ng and, thus, provide an empirical foundation for one of gravitation 
physics' fundamental components. Tests that reveal no violation of the Einstein equivalence 
principle force the form of the nongravitational Lagrangian density toward ones that ad- 
mit a metric representation, a representation in which a single metric tensor field couples 
universally |8j to all matter. 

The experimental tests of the Einstein equivalence principle that we analyze in this paper 
involve the propagation of polarized light through a background gravitational field. We are, 
therefore, concerned with only that part of the nongravitational Lagrangian density which 
governs the dynamics of the electromagnetic field and we base our analysis on a general 
model of this Lagrangian density introduced by Ni |J. One advantage of this approach is 
that the results of our analysis are easily specialized to yield the predictions of any theory 
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of gravity whose structure is encompassed by Ni's X9 formalism. The other advantage is 
that a broad consideration of conceivable nonmetric couplings provides a context in which 
to judge how different experimental tests of the Einstein equivalence principle compete with 
and complement one another and, so, to clarify the strengths and weaknesses of the empirical 
foundation of one of the most fundamental aspects of gravitation physics. 

The structure of the \g formalism is consistent with the basic tenets of the Dicke frame- 
work Q. Spacetime is a different iable manifold and matter fields and phenomenological 
gravitational potentials are tensor fields. Given the apparent smoothness of space and time 
and the view that physics should not depend in an essential way on a choice of coordi- 
nates, it is difficult to fault these axioms as a basis for phenomenology. Furthermore, it 
is assumed that as gravity is "turned off" the nongravitational Lagrangian density, Cng, 
reduces smoothly to the corresponding special relativistic Lagrangian density. The most 
general Lagrangian density consistent with these axioms and with assumptions of electro- 
magnetic gauge invariance, linearity of the electromagnetic field equations and the absence 
of couplings to derivatives of gravitational potentials is 

£em = -tt - X Q/37<5 F aj3 F lS , (3) 

1D7T 

where the electromagnetic field is related to a 4-potential in the usual way, F a p = Ap ;a —A a) p 
and where X aPl& is a tensor density that provides a phenomenological representation of 
the gravitational field. Note the analogy between ^ Q ^ <5 and a macroscopic electrodynamic 
constitutive tensor density. 

To see how the phenomenological field X aPl& is constructed from a specific theory's 
scalar, vector and tensor gravitational fields, the structure of the theory's electromagnetic 
Lagrangian density can be matched to that of Cem in Eq. (||). In metric theories of gravity 
it is constructed from the metric tensor alone, x a ^~ f& = \\/~ 9{g ai 9^ S ~ 9 a& 

Notice that x a ^ S inherits symmetries analogous to those of the Riemann tensor from 
the antisymmetry of F a p and from the symmetrical way in which it appears in Eq. ([$]). 
Consequently, x al3lS has twenty-one independent components. Ni says that his formalism 
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has twenty-one nonmetric degrees of freedom because, in general, one cannot exploit the 
formalism's covariance to transform away any of these independent component potentials. 
The freedom to redefine coordinates is exhausted by putting the formalism's g a p field into a 
standard form. This field represents the coupling between the gravitational field and scalar 
particles. 

III. GEOMETRIC OPTICS IN A BACKGROUND GRAVITATIONAL FIELD 

In this paper we consider the propagation of light through relatively weak gravitational 
fields that vary on length and time scales that are far longer than the light's wavelength and 
period respectively. We, therefore, adapt the coordinate system in which the Lagrangian 
density (||) is represented to the weak-field limit and we employ the methodology of geometric 
optics to analyze the light's propagation. 

In a Lorentzian coordinate system the electromagnetic Lagrangian density of special 
relativity has the form in Eq. (|3|) with -^ a ^ s = ^(j] ai r]^ s — rf-^rf 1 }. Consequently, in the 
weak field limit we can find quasi-Lorentzian coordinate systems in which the field x 01 ^ 5 
has the form 

x a^s = - fror/Tfe _ rj^rf-r) + 6x af3yS , (4) 

with 6x al3lS <C 1. Such a decomposition of covariant with respect to Poincare 

transformations provided that boost velocities too close to unity are not allowed. Note that 
this is a purely formal covariance. The background gravitational field may well single out a 
preferred frame ||. 

Geometric optics treats the propagation of locally plane electromagnetic waves. The 
amplitude and phase representation of such a wave, 

E = A E e i0 , B = A B e i0 . (5) 

is characterized by the fact that derivatives of the vector amplitudes Ae and Ab are small 
compared to derivatives of the rapidly varying phase function $. Since we are interested 
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in the propagation of high-frequency electromagnetic radiation through a weak background 
gravitational field that varies slowly in space and time, derivatives of the background field 
are also small compared to derivatives of $. We base the analysis that follows on the eikonal 
equation which determines local coordinate velocities of wave propagation. This equation 
is derived by inserting the representation (Q) into the electromagnetic field equations and 
ignoring all derivatives other than those of the phase function. Propagation equations for 
the vector amplitudes Ae and A# and equations governing post-geometric-optic corrections 
can be derived by taking systematic account of smaller derivatives but we shall not need 
these equations here. 

Neglecting derivatives of the slowly varying background field, the electromagnetic field 
equation that follows from Eq. (0) is 
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aM& F^ p = 0. (6) 



Defining electric and magnetic fields via F i0 = and = Sj^Bi and employing the 
decomposition (|]) of x a ^ 16 1 this can be written somewhat more transparently as 

V • E + terms proportional to 5% and E or B = 0, (7) 



and 



<9E 

V x B — - — h terms proportional to 5x and E or B = 0. (8) 
The electric and magnetic fields also satisfy 

VxE+ — =0, (9) 

and 

V-B = 0, (10) 

by virtue of the fact that F a p and, so, E and B are derived from potentials in the usual way. 

To derive the eikonal equation we insert the representation of a locally plane wave into 
the field equations (0) - (|10"D and neglect all derivatives other than those of the rapidly varying 



phase function $. Let denote the gradient of this function, k^ = <9 M $ = (d$/dt, V$) = 
(— cj,k), so that Eqs. (0) and ( JT0| ) become 

A B = ii^ (11) 

and 

k-A B = 0. (12) 

The latter of these implies that the magnetic field of a locally plane wave is transverse to 
the direction in which the wave propagates. Since Eq. (0) becomes 

k • Ag = terms proportional to 6x and or A#, (13) 

the electric field is guaranteed to be purely transverse only in the absence of gravity. To first 
order in the small <5\; Q ^ 7<5 that represent the weak background gravitational field, Eq. ([13]) 
expresses a wave's longitudinal electric field component in terms of its dominant transverse 
field components. To this same order, 0((5%), Eqs. (^T|) and ( |T3"D imply 

A E = —k x Ab + terms proportional to 6x and Ag, (14) 
k 



which expresses A^ in terms of the two independent components of A#, and Eqs. (§), ( |lT|) 
and fll~4]) imply the eikonal equation 

J 2 

(1 — — )A B = terms proportional to 5x and A B . (15) 
k 2 

Since the magnetic amplitude Ab has two independent components and since ou/ k is coor- 
dinate phase velocity, finding the two independent polarization states that propagate with 
well-defined phase velocities is a matter of solving a two-dimensional eigenvalue problem. 
Note that in the absence of gravity the right-hand side of Eq. (|15| ) vanishes implying, as one 
expects, u/k = 1 regardless of polarization. 

Before examining the explicit form of Eq. (^) it is convenient to make a 3+1 decompo- 
sition of 5x a ^ 1& into a set of SO (3) tensor objects. We define 
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f J = -5x mj i f j = -e jlm 5x° ilm , and ( ij = -e ilm e jpq 5x lmpq , (16) 

where e ijfc is the Levi-Civita antisymmetric symbol. This decomposition is covariant under 
rotations of the quasi-Lorentzian coordinate system in which Sx al3 ^ s is represented. It follows 
that the tensors Y J an d C 1 represent spatial anisotropy induced by the background 
gravitational field in a natural way. The symmetries of x al3lS and Sx al3lS imply that £ lJ and 
( l i are symmetric. Their twelve independent components combine with the nine of 7 U to 
account for the twenty-one nonmetric degrees of freedom of the X9 formalism. An analogous 
decomposition of x a ^ & yields tensors e u = 5^ + £ y , x y = ^ + J an d 7 U • 

Consider the propagation of light in the neighborhood of some event in spacetime. We 
exploit the covariance of our decomposition of 6x a/3 ^ s to rotate from the original quasi- 
Lorentzian (t, x, y, z) coordinate system in which the background gravitational field is rep- 
resented to a set of (t, x', y', z') coordinates in which the light propagates in the z' direction. 
In this new coordinate system A# has only x' and y' components so Eq. ( |i5| ) reduces to the 
system of two equations, 

{1- ( £)A$ = AA%-BA%, (17) 

{1-£)A% = -BA*+CA%. (18) 

The coefficients A, B and C depend on location in spacetime and on the direction in which 
the wave propagates, the z' direction. This dependence is implicit in the expressions for A, 
B and C in terms of the values at the event in question of the tensor components , C*'- 7 ' 
and Y' 3 ' m t ne {t,x',y',z') coordinate system, 

A=^ -2-f 1 ' -C 1 ' 1 ', (19) 
i3 = ^' 2 ' + ( 7 2 ' 2 '- 7 1 ' 1 ') + C 1 ' 2 ', (20) 

and 
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C = C 1 ' 1 ' + 27 1 ' 2 ' -C 2 ' 2 '. (21) 

The matrix defining the structure of the right-hand sides of Eqs. (|17|) and ([18]) is real- valued 
and symmetric when (5^ a ^ 7<5 is real. Its eigenvalues are 

\± = ^±y(A-Cy + 4B*. (22) 

The corresponding eigenvectors define the polarization states that propagate with well- 
defined coordinate phase velocities c± = 1 — |A± + 0(Sx 2 )- 

The effects of a nonmetric background field on physical standards of length and time are 
such that a local observer who uses rulers and clocks to measure velocities of light having 
the polarizations singled out by the field will generally obtain results that differ from the 
coordinate velocities c + and c_ at 0(5x)- To that order, however, the fractional difference 
between the velocities that an observer at rest in the (t, x', y', z') coordinate system does 
measure will be the same as the fractional difference between c + and c_. We let 5c/ ' c denote 
the magnitude of this dimensionless local observable, 

^ = y(A-Cy + ^ + 0(5x 2 ). (23) 

A practical way to measure or to impose an upper limit on this magnitude is to search for 
the effect such birefringence has on the polarization of light that propagates through the 
background gravitational field. 

We use Stokes parameters based on the linear polarization states singled out by the 
background field to characterize the polarization of light that propagates through a local 
observer's laboratory. The intensity / and the degrees of polarization Q/I, U/I and V/I 
provide a particularly convenient representation of the light's state. We may think of the 
light as an incoherent superposition of appropriate intensities of unpolarized light, linearly 
polarized light with Q/I = ±1 and elliptically polarized light with (U 2 + V 2 ) 1 / 2 /I = ±1. 
Propagation through the observer's laboratory only affects the state of the last of these 
components. 
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The unpolarized component is itself an incoherent superposition of equal amplitudes of 
the basis polarization states while the linearly polarized component consists purely of one 
or the other of the basis polarization states. The shift in the relative phase of waves in the 
basis states that accumulates as these waves propagate with different phase velocity has no 
effect on the polarization of either of these components. On the other hand, the elliptically 
polarized component is a coherent superposition of equal amplitudes of the basis states so 
its polarization is sensitive to the relative phase of these waves. Indeed, the relative phase 
determines this component's U/I and V/I degrees of polarization, V/U = tan($ + — $_). As 
light propagates during a coordinate time interval dt this relative phase shifts by uj(5c/c)dt, 
where u is the light's angular frequency. It is this shift that can affect light's polarization. 

We piece together the cumulative effect of propagation through a background gravita- 
tional field on a signal's polarization from such local effects. This is particularly easy to 
do when the orientation of the basis polarization states singled out by the background field 
is fixed along light rays, as is the case for light propagating through the essentially static, 
spherically symmetric background fields that we consider in the next section. A ray thread- 
ing through a spherically symmetric field lies in a plane and, for the fields we consider, there 
is a fixed orientation between this plane and the A# amplitudes of the polarization states 
singled out by the fields. In such cases the relative phase shift that accumulates between 
waves in the basis states is simply 

A<S> = u J^dt + 0(5x 2 ), (24) 

where, since the local value of 5c/ c is explicitly 0(Sx), the integration runs along the un- 
perturbed ray that runs straight through the quasi-Lorentzian coordinate system with unit 
speed from source to receiver. The integrand in Eq. (pi]) is the local value of 5c/ c for the 
ray's direction. 

One consequence of the relative phase shift is that a signal whose elliptically polar- 
ized component has V/U — when emitted will have some degree of circular polarization 
when it is received. The degree of circular polarization received is proportional to a small 
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accumulated shift, 

V/U « A$. (25) 

We base our interpretation of galaxy polarization data in the next section on this expression. 

We have seen that a nonmetric background field induces birefringence if and only if 
(A — C) or B in Eq. (^) is nonzero for some direction of wave propagation at some event 
in spacetime. Before we turn to the derivation of empirical limits on such birefringence let 
us consider how such limits constrain the form of the coupling between gravitational and 
electromagnetic fields. 

Note that the expressions (A — C) and B involve only the symmetric, trace-free parts 
of the tensors £ lJ , and 7 * J '. This is most easily appreciated by writing (^4 — C) and 
B in terms of the spherical components of these tensors. In terms of components in the 
(t, x', y', z') coordinate system they are 

a-c = A [(e (?) + e g, } + 2<( 7 c?) - 7 g,) + (ci 2) + cS)] (26) 

and 

B = ~ &) + 2(# + 7 i|) + - Ci|)]. (27) 

Only / = 2 components appear. Expressions for (^4 — C) and B in terms of components in 
the original (t, x, y, z) coordinate system follow from the transformation law for spherical 
tensor components |10[ , for example, 

where 0, ^ and i/j are the Euler angles specifying the rotation from (t, x, y, z) to (t, x', y', z'). 
From the form of this transformation law and of Eqs. (p6|) and (^7]) we conclude that the 
absence of gravity- induced birefringence implies + = and 7 ^ = and, so, that 
observations imposing limits on the strength of gravity-induced birefringence constrain ten 
of the X9 formalism's twenty-one nonmetric degrees of freedom. This conclusion and the 
expression ( |2"B| ) for Sc/c are in accord with remarks of Ni [0 for which the preceding analysis 
provides a detailed justification. 
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IV. GRAVITY-INDUCED BIREFRINGENCE AND THE POLARIZATION OF 

GALAXIES 



The preceding analysis examined nonmetric couplings between the gravitational and 
electromagnetic fields and established which single out pairs of linear polarization states 
that propagate through a gravitational field with different phase velocities. Using Stokes 
parameters based on those states, we found that such gravity-induced birefringence causes 
[/-polarized radiation to develop a degree of circular polarization as it propagates. 

Here, we interpret polarization measurements of extragalactic sources cataloged by Weiler 



and De Pater |LT] in light of this prediction and derive a precise new constraint on the 
strength of birefringence that could be induced by the Galaxy's gravitational field. Note 
that even a cursory examination of the data indicates that such birefringence cannot be 
strong. If it were, one would expect to measure substantial circular and linear polarization 
of radiation from extragalactic sources with about the same likelihood. This is not the case. 
Measurements reveal significant circular polarization of radiation from extragalactic sources 
only rarely while they frequently reveal substantial linear polarization. 

As a first approximation, we treat the Galaxy's gravitational field as static and spherically 
symmetric. This is not unreasonable, especially if one accepts the existence of a substantial 
dark-matter halo. Note that the assumed spherical symmetry does not preclude the detection 
of gravity- induced birefringence. Radiation reaching the Earth from an extragalactic source 
generally does not propagate radially through the Galaxy's field because the Earth lies 10 
kiloparsecs from the Galaxy's center. A ray along which such radiation travels lies in a plane 
defined by the locations of the source, the Earth and the center of the Galaxy. 

At each point along a ray the only direction that the Galaxy's spherical field can single 
out is the radial one. We, therefore, find it useful to introduce at each of these points a local 
quasi-Lorentzian (t, x, y, z) coordinate system oriented so that its z axis is radial and its x 
axis lies in the ray's plane. We do so because the representation of the spherical tensors 
introduced in the preceding section is simple in these local systems. Specifically, 
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and ■yffl are nonzero only for m — 0. The nonzero components depend only on the radial 
coordinate of the point along the ray. 

At each of these points we can rotate the local (t, x, y, z) coordinate system about its 
y axis through an angle 9 to obtain a local (t,x',y',z') coordinate system in which the ray 
runs in the z' direction. Equation (|23|) expresses the local value of 5c/ 'c for points on the 
ray in terms of (A — C) and B which are, in turn, expressed in terms of the local (t, x', y', z') 
components £±l>, C±2' an d 7±y by Eqs. fl26|) and (p7|). Since the Euler angles of the rotation 
from (t, x, y, z) to (t, x', y', z') are = 0,6' and ip = 0, the transformation law ( p8| ) implies 

£l = sm\6)tf\ (29) 

with the same relationship holding between (±1, and ffl and between and 70 2 "*. For 
any spherically symmetric gravitational field, we conclude that 5c/ c is proportional to the 
square of the sine of the angle by which a ray's direction differs from the radial one. 

The orientations of the A# amplitudes of the polarization states singled out along a ray 

(2) (2) (2) 

are determined by the relative magnitudes of the (t,x,y,z) components Q , Q an d 7o • 
The form of Eqs. ( p6| ) and ( p?D and the transformation law ( p9"D imply that B is proportional 
to 7q 2 ' ) while {A — C) is proportional to + Co 2 '')- I n the event that (A-C)j^O and B = 0, 
Eqs. ( |17D and (|l^) imply that for the polarization states singled out by the Galaxy's field 
A# will lie in and perpendicular to the ray's plane. On the other hand, when (A — C) = 
and B 7^ 0, Eqs. (17) and (18) imply the Ab singled out will be oriented at 45 degrees 
to either side of the normal to the ray's plane. Cases like these, in which the polarization 
states singled out by a background field are fixed along rays, are special. Generally, their 
orientations will change along a ray as the relative magnitudes of (q 2 ^ and % vary. 

The effect that propagation through the Galaxy's gravitational field has on the polar- 
ization of light from extragalactic sources clearly depends on the global structure of that 
gravitational field. We choose to compute the effect for an NGT model of the Galaxy's field 
studied by Krisher || so that we can compare the precision of the constraint we impose 
on gravity- induced birefringence with the one he imposed. This NGT field has £q — 0, 
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(q 2 ^ = L 2 (r) and Jq = 0, where the function L 2 (r) specifies the structure of the antisym- 

(2) 

metric tensor part of the Galaxy's nonsymmetric field. The fact that only Q is nonzero 
is a consequence of the nonmetric coupling between the nonsymmetric field and the elec- 
tromagnetic field having been tuned to satisfy constraints imposed by other tests of the 

(2) 

Einstein equivalence principle 12|. Since only Q is nonzero, the orientations of the basis 
polarization states singled out by the Galactic field are fixed along rays with their mag- 



netic amplitudes lying in and perpendicular to ray planes. Consequently, Eq. (|24]) gives 
the relative phase shift that accumulates between waves in these polarization states as they 
propagate. We have 5c/ c = |L 2 sin 2 (6 l ), which is consistent with the results of Gabriel |]|| 
et al. since they denote L 2 by Q. 

In Krisher's model the antisymmetric part of the Galaxy's nonsymmetric gravitational 
field is generated by a uniform density, spherical halo of weakly interacting particles. Its 
radius, R, is 25 kiloparsecs. The function L 2 (r) is proportional to the halo's net NGT 
charge, I 2 . We have L 2 = l 4 /r 4 for radii greater than R and L 2 = l 4 r 2 /R 6 for smaller 
radii. The NGT charge, I 2 , can be expressed in terms of a coupling constant f 2 and the 
density of halo particles, n c , which Krisher takes to be 0.1 per cubic centimeter. Specifically, 
I 2 = (4nR 3 /3)n c f 2 . 

The geometry of a ray that runs from an extragalactic source to the Earth is shown in Fig. 
1. The spherical symmetry of the Galaxy's field implies an axial symmetry of the relative 
phase shift (M)- It depends only on the angle j3. The 10 kiloparsec distance between the 
Earth and the Galaxy's center is denoted by d. The ray's impact parameter is b = dsin(j3). 
To compute the phase shift (24), we parametrize the rays so that the point r = b, which 
may lie on an extension of the ray, corresponds to t — and we suppose that the ray 
starts from r = 00 at t = — 00. It follows that the ray encounters the halo's surface at 
t = t R = -(R 2 - & 2 ) 1/2 and the Earth at t = t E = -dcos(/3). Note that r 2 = (b 2 + 1 2 ) and 
sin 2 (6>) = b 2 /(b 2 + t 2 ) along the ray. 

From these facts and the expressions for L 2 inside and outside the halo we conclude that 
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A$ - = — jL (y^F = + 2ft* (*» + 1*) + 2y arctan ^ ))L - (30) 

and 

ul 4 b 2 ft* , ul 4 b 2 , 

Equation (|30D corrects a misprint in Krisher's expression for AQ ext . More seriously, Krisher's 
numerical estimate of A$ ex4 indicates that rather than £r was used as an upper limit of 
integration. This contribution to A$ was, therefore, significantly overestimated. In Fig. 2 
we plot conveniently normalized functions F ext {(3) and F int ((3) defined by A$j = (ul 4 /d 3 )Fi. 

Data from the catalog of Weiler and De Pater impose a sharp new limit on the amplitude 
of A$(/?) and, so, on the magnitudes of I 2 and f 2 . 

From the position of a source in the sky and a bit of geometry we obtain both the angle (3 
for the source and the orientation of the polarization states singled out for it by the Galaxy's 
field. We use this latter information to express the source's cataloged polarization in terms 
of Stokes parameters based on the polarization states singled out. Since the relative phase 
shift A$ converts U polarized light to V polarized light, we calculate V/U for sources with 
accurately measured U / 1 and V/ 1 and for which the U / 1 degree of polarization is substantial, 
accounting for more than half of the observed linear polarization. The substantial U/I degree 
of polarization means that a given relative phase shift will induce a relatively large V j I and 
it also assures that V/U will be a well-behaved indicator of A$ via Eq. (|25|) where we make 
the conservative assumption that a nonzero average of V/ U for sources in an interval around 
(3 represents the effect of the phase shift A<3>(/3). Figure 3 is a plot of the V/U as a function 
of P for data taken at 5 GHz, the highest frequency included in the Weiler-De Pater catalog. 
The average of V/U for these sources is consistent with zero. We estimate limits imposed 
by this data by noting that 0.12 is the sample standard deviation of V/U for the 5GHz data 
having f3 > 60 degrees. Equation (|25|) imposes this as a limit on the mean of A$ for (3 in 
this range. This limit and the expression in terms of f 2 for the mean of A$(/3) in the range 
from 60 to 90 degrees implies f 2 < 1.4 x 10~ 34 cm 2 . When the means and sample standard 
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deviations of V/U for lower frequency data are accounted for this constraint is sharpened 
slightly, 

f c < 1.3 x l(T 34 cm 2 . (32) 

This limit is more than 10 4 times tighter than the one Krisher derives, once his overesti- 
mate of the exterior contribution to A<3> is corrected for. This increase in precision is easily 
understood. Krisher argued that NGT-induced birefringence would affect the polarization 
structure of pulses from pulsars in the Large and Small Magellenic Clouds and reckoned that 
an effect would have been noticed if the arrival times of pulse components in the polarization 
states singled out by the Galaxy's field were shifted relative to each other by as much as a 
millisecond. In our case, we use a kind of interferometry to discern a much smaller relative 
time delay. The highest frequency data that we analyze imposes a limit of 0.12 radians on 
the relative phase shift between the polarization states singled out by the Galaxy's field. 
We are, therefore, resolving a fraction 0.12/27T of a cycle at the 5 GHz observation fre- 
quency. This corresponds to a relative time delay of 4 x 10~ 12 seconds. Our time resolution 
is, therefore, more than 10 8 times Krisher's. Since, as Krisher remarks, constraints on the 
magnitudes of I 2 or f 2 are proportional to the square root of the time interval resolved, the 
precision of our constraint ([32|) is accounted for. 

V. SUMMARY AND CONCLUSIONS 

The local isotropy of space is one consequence of the Einstein equivalence principle. 
Atomic physics experiments of Hughes-Drever || type are the classic tests of this predic- 
tion. In this paper we have shown that searching for effects of propagation through a 
gravitational field on the polarization of light provides a precise, qualitatively different test 
for spatial anisotropy, specifically spatial anisotropy that induces a birefringence of space. 
In the context of the X9 formalism such anisotropy is associated with ten of twenty-one 
nonmetric degrees of freedom. 
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The significance of any test of the Einstein equivalence principle lies in the constraint it 
imposes on the form of the coupling between gravitational and matter fields. Constraints on 
gravity- induced birefringence force the structure of the electromagnetic Lagrangian density 



That the phenomenological symmetric tensor potential g a ^ and scalar potential can rep- 
resent the X9 formalism's remaining eleven nonmetric degrees of freedom is clear upon 
counting components. The fact that this Lagrangian density predicts no birefringence in 
the geometric optics limit is clear from the metric form of the coupling between g a @ and 
the electromagnetic field and from the identification of as the 7q°' ) part of 5x a ^ 1& '■ Note, 
however, that Carroll and Field ]13[ have shown that at post-geometric-optics order the 
coupling does induce a dispersive optical activity which they have constrained using galaxy 
polarization data. 

Two features are primarily responsible for the precision of Einstein equivalence principle 
tests based on the interpretation of polarization data for extragalactic sources. First, they 
focus on global consequences of Einstein equivalence principle violations. They constrain 
shifts in relative phase that accumulate from minute local effects rather than seeking the 
local effects directly. Second, they are interferometric in character. Their sensitivity to 
accumulated phase shifts comes from exploiting the coherence inherent in polarized light. 
Note that the sharp new constraint on NGT derived above is far beyond the reach of local 
Hughes-Drever experiments, as indicated by an estimate analogous to one that Gabriel |J 
et al. make of local anisotropy induced by the Sun's field in NGT. 

Sharper tests of the Einstein equivalence principle are all but certain to result from 
surveys of extragalactic sources that yield linear and circular polarization measurements at 
frequencies higher than the 5 GHz limit of the Weiler-De Pater catalog. 
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FIGURES 

FIG. 1. Points G, E, and S depict the locations of the Galaxy's center, the Earth and an 
extragalactic source. A halo of radius R = 25kpc surrounds the galaxy. The Earth lies at a 
distance d = Wkpc from G. Light from the source S propagates to the Earth along the ray SE. 
The relative phase shift that accumulates depends only on the angle (3 or, equivalently, on the 
impact parameter b = dsin(/3). At each point P of SE we denote the angle between SE and the 
radial line GP by 9{P). 

FIG. 2. The function of (3 depicted by the solid curve is proportional to the relative phase 
delay. It is the sum of the internal and external contributions, Fi nt ([3) and F ext ([3), defined in the 
text and denoted here by dotted and dashed curves respectively. 

FIG. 3. Plot of V/U versus (3 for 5GHz data from the Weiler-De Pater catalog [11]. The 
apparent clustering of data is a consequence of including data taken at different times for some 
sources. 
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